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Abstract 

It is rigorously shown that the two-dimensional Hubbard model with nar- 
row bands (including next nearest-neighbor hopping, etc.) does not exhibit 
d^2_y2-waye pairing long-range order at any nonzero temperature. This kind 
of pairing long-range order will also be excluded at zero temperature if an 
excited energy gap opens in the charge excitation spectrum of the system. 
These results hold true for both repulsive and attractive Coulomb interac- 
tions and for any electron fillings, and are consistent with quantum Monto 
Carlo calculations. 

PACS numbers: 74.20.-z, 05.30.-d 
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Considerable experimental evidence shows that the dominant symmetry of the pairing 
order parameter in high temperature superconductors may be dx'i_y2-waNe (see, e.g., Ref. 
|I|] for a review). Though some disputes regarding this challenging issue still remain, a 
great number of people arrive at such a consensus that exploring the possibility of d-wave 
superconductivity in strongly correlated electrons would be quite useful towards ultimately 
successfully explaining high temperature superconductivity, thus resulting in numerous stud- 
ies on this subject. As a matter of fact, owing to its apparent simplicity the two-dimensional 
(2D) Hubbard model naturally becomes an actively debating focus in recent years, as it is 
widely thought to provide a simple model to interpret some essential features relevant to 
the physical properties of Cu02 planes in the cuprate oxides. In spite of intense efforts 
being made both numerically and analytically 0,^, however, a basic question whether or 
not the dx'i-y-i-waNe pairing long-range order (LRO) in the 2D Hubbard model exists, is still 
inconclusive. Actually, at energy scales and lattice sizes accessible to numerical simulations 
no definite sign of d^^^yi superconductivity has been detected in this model, while some 
analytical works using different approximations appear to suggest positive answers 
thereby leaving some controversies and ambiguities to be resolved. To clarify them, rigorous 
results are particularly needed at this stage. 

In this paper, based on Bogoliubov's inequality we show rigorously that the 2D Hubbard 
model with narrow bands (including next nearest-neighbor hopping, etc.) does not exhibit 
dx2_y2-'wave pairing LRO at any nonzero temperature. This kind of pairing LRO will also be 
excluded if a gap opens in the charge excitation spectrum of the system. These results hold 
true for both repulsive and attractive Coulomb interactions and for any electron fillings. 
Combining with other known exact results, one would conclude that the 2D Hubbard model 
might not have the right stuff for describing superconductivity in the cuprate oxides in this 
sense provided that the superconducting mechanisms in these materials are supposed to be 
due to condensation of either s-wave Cooper, or generalized rj or dx^-y^-wave electron pairs. 
The present observations are consistent with quantum Monto Carlo results. 

Let us start with some preliminary definitions. The dx^-y^-wave pairing operator (like 
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the Cooper pairing operator) can be defined as 



= II(C0S h - COS fcy)4^clki = IY1 fi^HA+Sl^ = (^d )^ (1) 

k ^ r,S 

where we have used the Fourier transform of an electron operator Cko- 

Cka = exp(ik ■ Kj)cj^, Cj„ = exp(-ik ■ Rj)ck^, (2) 

and f{d) = +1 (—1) for S = ia^. (ia^), and zero otherwise, where a^,, ay are unit vectors 
connecting nearest-neighbor sites, M is the number of lattice sites, Rj (r) is the position 
vector of jth (rth) site, and a denotes spin. In the following we for simplicity take S = ia^,, 
±ay. According to Bogoliubov [0, when one studies a degenerate state of the statistical 
equilibrium, one should first remove the degeneracy by introducing a symmetry-breaking 
field, and then turn to investigate the so-called quasi-averages involved. On account of this 
reason we define the dx2-y2-wa:Ve pairing order parameter per site as 

A+ 

g= lim lim (-f), (3) 

where (■ ■ ■) stands for the thermal average in a grand canonical ensemble, and u is the 
amplitude of a U(l) symmetry-breaking field. Note that the two limit processes are non- 
interchangeable. The nonvanishing of g, namely g ^ means the existence of dx2-y2-waye 
pairing LRO, while g = gives the converse result. 

We now consider a general Hubbard model with narrow bands on a periodic lattice in 
the presence of a U(l) symmetry-breaking field. The Hamiltonian reads 

i,j c i i 

-u{Aj + A^), (4) 

where the sum on i and j can run over all M lattice sites, (qo-) is the creation (an- 
nihilation) operator for an electron at site i with spin a, the on-site Coulomb interac- 
tion Ui and the chemical potential fii are allowed to be position-dependent for generality, 
and rijo- = cl^Ci^j, the number operator of electrons. T(Rj — R^), which has a property 



T(R.j — Hj) = T*(Rj — Hi), is the local overlap integral which designates the energy bands 
of the model. In fact we require that T(Ri — Rj) survives only for short-r anged overlapping 
in the present case. The last term in Eq.(^) introduces the effect of the U(l) symmetry- 
breaking field, where u is an infinitesimal quantity, and are given by Eq.(|l]). 

As we attempt to explore the possibility of dx2-y2 pairing LRO in Eq.(^), in the following 
we shall use Bogoliubov's inequality 

\{[C\A^])?<^{{A,A^]){[[C,HIC^]), (5) 

for any quantum-mechanical operators A and C, and with the inverse temperature (3 = 1/T 
{ks = 1), where [, ] and {, } are the usual commutator and anticommutator, respectively. 
This inequality was used to exclude the possibility of magnetic LRO in the Heisenberg ^ 
and Hubbard |llO| , pl]| models as well as superfluidity in Fermi liquids in one and two 



dimensions at nonzero temperature. Quite recently, this inequality was also used to exclude 
the possibility of s-wave Cooper pairing and generalized rj pairing LRO in the ID and 2D 
Hubbard models with narrow bands at nonzero temperature It should be mentioned 
that the SU(2) Lie algebra obeyed by the relevant operators (e.g., spin operators. Cooper 
pairing operators and r] pairing operators, etc.) plays a key role in applying this inequality 
to the above-mentioned cases. However, one may observe that the d^^-yi-waNe pairing 
operators defined in Eq.(|l]) do not obey the SU(2) symmetry, which makes Bogoliubov's 
inequahty not directly applicable to the present stated in Ref. |13|. Fortunately, 

this difficulty can be overcome by noting the simple fact that /(<5) in Eq.(|l|) takes values 
either +1 and —1 or zero so that we can decompose the dx2_y2 pairing order parameter per 
site, g, into four terms each of which obeys the SU(2) algebra. It is this property that makes 
it possible for applying inequality (||) to our case |T^. We would like to mention here that 
although one can find the standard derivations in Refs. p|- p!l| , p!3| , for reader's convenience 



and for this paper being self-contained, we shall below intend to present our calculations in 
some detail. 

We define the following operators 



Vt = clA+cV = Cr+«iCrT, Vr = ^(^rj + «r+«i - 1), (6) 

with an arbitrary constant vector a on the lattice. It can be verified that they satisfy 
The Fourier transforms of fj operators, hke spin operators in Refs. ||9|JTT[], are defined by 



vf'' = 7^ E exp(-ik ■ r)r/±'^(k), f/±'"(k) = ^ exp(ik ■ r)f/±'^ (8) 

k r 

which comply 

[r7+(k),77-(k')] =2r7^(k + k'), [r?±(k), ^^^(k')] = T^±(k + k'). (9) 

Eqs-dl) come from Eqs.(^ and (|). 

With these definitions we choose A = ?7+(— k— Q) and C = f/^(k) in for our purposes. 
After some algebra for the double-commutator one gets the inequality, 

■ik),H],fi^{-k)]) < i5:|T(R,)||cos(k-R,)-l||^e^'^'-^'(cLcuv)| 



TV 

<jT.\Tm\Rlk' + 2\u\.\{Aj)\. (10) 

i 

In the derivation of this inequality we have used the property of translation invariance and 
such a few simple facts as J^k'cric'i'a'^ii'a) = N, the total number of electrons, 1 — cosx < 
and (A^) = (A^)''". Substituting inequality ([T0| ) into (H) we have 

(-Q - riQ H- k)}> > , Jj^tff(A" /M)i) ■ (") 

where J^.,Af(Q,a) = (f/+(Q))/M, and ^ = {N/M)E^\T{'Ri)\Rj /4. Since the r(Ri)'s 
are the matrix elements of the overlap integral between Wannier functions which de- 
crease rapidly with distance for strongly correlated electrons, the summation J2i \T(Ri)\Il'^ 
is well defined. For the single-band Hubbard model as well as one with next nearest- 
neighbor hopping integral the values of J2i\T(Ri)\Ilf on a hypercubic lattice can be 



found in Ref. |T^. Summing both sides of the above inequahty over k and noting that 
{1/M)Y.MAA*}) = EAl¥,rV}} = Er{|l - (n,t - nr+„i)'l> < M, we obtain 

Now we take the thermodynamic hmit, i.e., N ^ oo and M — ^ oo with the ratio N/M 
fixed. Then the sum on k in ([T2|) can be replaced by the integral over the first Brillouin 
zone. Suppose that fco is the distance of the nearest Bragg plane from the origin in k space. 
Then we obtain for small z/ the following inequalities 

|^.,oo(Q,a)|<(e/3')l/>^/^ (ID) (13) 

Ud 1 

|^.,oo(Q,a)|< (2D) (14) 

V TT I m |z/| |-^/^ 

where we have used limjv/^oo I /^) I < 2 in (|1^). Inequalities (|13D and (|14D tell 
us that |jFo oo(Q, Q:)| = for any nonzero temperature as the U(l) symmetry-breaking 
field is turned off (z/ 0+). (Recall that the thermodynamic limit has been taken 
before we remove the U(l) symmetry-breaking field.) This means that ^o,oo(Q)Q;) := 
limj,^o+ liinAf^oo((l/^^) J2r exp {iQ ■ r)4-^4_^„|)j,^M = for any possible Q and a for T > 
in the ID and 2D Hubbard models defined in Eq.(^. On the other hand, we note that the 
dx2-y-2-vfave pairing order parameter per site can be rewritten as g = (1/2) J2s /(^)-^o,oo(0, d) 
because the existence of J^o,oo(0, S) has been proved, where we have set Q = and ct = S. 
Here use has been made of the well-known theorem: lima ± lim6 = lim(a ± b) if lima and 
lim6 exist. Since ^o,oo(0, <5) = for T > and f{S), by virtue of its definition, takes either 
±1 or zero, we finally have g = ioi T > 0. Consequently, we have proved that the 2D 
Hubbard model with narrow bands, defined by Eq.(^, does not show dx2-y2-wave pairing 
LRO at any nonzero temperature. 

When P/2 is replaced by 1/Egap in inequality (H), where Egap (> 0) is an excitation 
energy gap between the lowest excited state and the ground state of the system, still 
holds true at zero temperature. We refer to Refs. [|I^,|r^ for detail discussions. Therefore, all 
the above analyses can apply to the case at zero temperature except P/2 replaced by 1/Egap. 
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Then we could conclude that if an energy gap opens in the charge excitation spectrum of 
the system (|), there will also be no dr^2_y2-wme pairing LRO in two dimensions at zero 
temperature. 

A few remarks are in order, (i) The derivations above presented are valid for both 
repulsive and attractive on-site Coulomb interactions and for any electron filling fraction, 
(ii) The method used in this paper can be readily extended to exclude the possibility of 
extended s- wave [with pairing operator A^_^ = I]k(cos /c^^ + cos A;y)cj^|cl,^J superconductivity 
in the ID and 2D Hubbard models at finite temperatures. In the Hilbert subspace without 
doubly-occupied sites this method might be applied to the t-J model as well. It should 
be emphazied that the method works only for low dimensions (ID and 2D), not for three 
dimensions, (iii) The conclusions drawn in this paper are quite consistent with the results 
from numerical simulations (e.g., quantum Monte Carlo calculations) in the 2D Hubbard 
model. We refer to Refs. for excellent reviews. The practical situation is that although 
most numerical works show some tendencies favouring d^2_y2 superconductivity, but no 
definite sign of LRO has been detected, as pointed out in Refs. |]TB|,|T7|. We would like to 
mention here that the present conclusions are not incompatible with the quantum Monte 
Carlo observations that the long-tailed enhancements in the dx2-y2 pairing correlation near 



half-filling or the exhibition of dx2-y2-\ike pairing fluctuations at low temperatures |T8 
are detected in the 2D Hubbard model. Recent analytic and quantum Monte Carlo results 
also show that the 2D Hubbard model with next nearest-neighbor hopping integral does not 
exhibit any definite sign of s-wave and d-wave superconductivity ||T9|-p^, consistent with 
the present exact result, (iv) The nonexistence of s-wave Cooper pairing and generalized r) 
pairing LRO at finite temperatures in the ID and 2D Hubbard models with narrow bands 
has been proved in Ref. Combining these exact results one would conclude that 

the 2D Hubbard model might not have the right stuff for explaining high temperature 
superconductivity in the layered cuprate oxides if the superconducting mechanisms in these 
materials are supposed to be due to condensation of one of the above-mentioned electron 
pairs, as a successful t heory should describe unifyingly the properties not only at zero 



temperature but also at finite temperatures, (v) To choose a proper model which could 
exhibit c?a.2_y2-wave pairing LRO in 2D, one may consider those which contain electron 
interactions not commuting with the operator X]r6xp(zk ■ r)?7,r, like one investigated in Ref. 
p!6| , because in this way Bogoliubov's inequality becomes ineffective. Another possibility is 
that the coupling between layers in the cuprate oxides might ought to be considered, which 
could also make Bogoliubov's inequality ineffective. 

To summarize, we show rigorously that, by means of Bogoliubov's inequality, the 2D 
Hubbard model with narrow bands (including next nearest-neighbor hopping, etc.) does 
not exhibit (i^2_y2 wave pairing LRO at any nonzero temperature. This kind of pairing LRO 
will also be excluded if an excited energy gap opens in the charge excitation spectrum of 
the system. These results are valid for both repulsive and attractive Coulomb interactions 
and for any electron fillings. Combining with known exact results obtained previously one 
would conclude that the 2D Hubbard model might not have enough right stuff for describing 
superconductivity in the cuprate oxides if the superconducting mechanisms in these materials 
are supposed to be due to condensation of one of the aforementioned electron pairs. The 
present observations are consistent with quantum Monto Carlo results. 

ACKNOWLEDGMENTS 

One of authors (GS) is grateful to the Department of Applied Physics, Science University 
of Tokyo, for the warm hospitality, and to the Nishina foundation for support. This work 
has been supported by the CREST (Core Research for Evolutional Science and Technology) 
of the Japan Science and Technology Corporation (JST). 



8 



REFERENCES 



*0n leave from Graduate School, Chinese Academy of Sciences, Beijing, China. Elec- 
tronic address: gsu@ap.kagu.sut.ac.jp 
'''Electronic address: msuzuki@ap.kagu.sut.ac.jp 

[1] Z.-X. Shen and D.S. Dessau, Phys. Rep. 253, 1 (1995) and references therein. 

[2] P.W. Anderson, D. Pines, D.J. Scalapino and Z.X. Shen, Phys. Today 47, No.2, 11 
(1994). 

[3] D.J. Scalapino, Phys. Rep. 250, 329 (1995); D. Van Harlingen, Rev. Mod. Phys. 67, 
515 (1995); and references therein. 

[4] E. Dagotto, Rev. Mod. Phys. 66, 763 (1994) and references therein. 

[5] D.J. Scalapino, in Does the Hubbard model have the right stuff?, edited by R.A. Broglia 
and J.R. Schrieffer, Proceedings of the International School of Physics, Enrico Fermi, 
Course CXXI (North- Holland, Amsterdam, 1994). 

[6] Note that without loss of generality we adopt this simple form for the d^i-yi pairing 
operator (see, e.g. Ref. @]). If a more symmetrical form, like one exploited in Refs. [^||, 
is used, the present method also works, and the conclusion drawn in this paper still 
remains true. 

[7] N.N. Bogoliubov, Physica 26, SI (1960). 

[8] N.N. Bogohubov, Phys. Abh. Sowjetunion 6, 1, 113, 229 (1962); see also H. Wagner, Z. 
Phys. 195, 273 (1966). 

[9] N.D. Mermin and H. Wagner, Phys. Rev. Lett. 17, 1133 (1966). 

[10] M.B. Walker and Th.W. Ruijgrok, Phys. Rev. 171, 513(1968). 

[11] D.K. Ghosh, Phys. Rev. Lett. 27, 1584 (1971); G.S. Uhrig, Phys. Rev. B 45, 4738(1992). 



9 



[12] P.C. Hohenberg, Phys. Rev. 158, 383 (1967). 

[13] G. Su, A. Schadschneider and J. Zittartz, Phys. Lett. A230, 99 (1997). 
[14] This simple fact was unfortunately overlooked in Ref. [|T^|. 



[15] Assa Auerbach, Interacting Electrons and Quantum Magnetism ( Springer- Verlag, New 
York, 1994). 

[16] F.F. Assaad, M. Imada and D.J. Scalapino, Phys. Rev. Lett. 77, 4592 (1996). 
[17] K. Kuroki and H. Aoki, preprint |cond-mat /9 70 72 77 . 



[18] D.J. Scalapino, preprints |cond-mat / 960621^ , 9606218 



[19] R. Lacaze, A. Morel, B. Petersson and J. Schroper, preprint |cond-mat/96 10193 



[20] A.F. Veilleux, A. Dare, L. Chen, Y. Vilk and A. Tremblay, Phys. Rev. B 52, 16255 
(1995). 

[21] T. Husslein, I. Morgenstern, D.M. Newns, P.C. Pattnaik, J.M. Singer, and H.G. Matut- 
tis, Phys. Rev. B 54, 16179 (1996). 

[22] S. Zhang, J. Carlson and J.E. Gubernatis, Phye. Rev. Lett. 78, 4486(1997); Phys. Rev. 
B55, 7464 (1997). 

[23] T. Koma and H. Tasaki, Phys. Rev. Lett. 68, 2348 (1992). It can be noted that the 
method presented in this paper could be extended to cover the d-wave pairing case if 
a similar argument to ours (i.e., dividing the d-wave pairing operators into four terms) 
is taken into account. Our presentation, however, is simpler and accessible to general 
reader. 



10 



